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The semiclassical limit for Dirac particles interacting with a static gravitational field is investi- 
gated. A Foldy-Wouthuysen transformation which diagonalizes at the semiclassical order the Dirac 
equation for an arbitrary static spacetime metric is realized. In this representation the Hamiltonian 
provides for a coupling between spin and gravity through the torsion of the gravitational field. In 
the specific case of a symmetric gravitational field we retrieve the Hamiltonian previously found by 
other authors. But our formalism provides for another effect, namely, the spin hall effect, which was 
not predicted before in this context. 

O ' PACS numbers: 

o: 

;-( ' INTRODUCTION 

'^^ . The dynamics of quantum particles in a gravitational field is an important topic, at the crossroad of particle Physics 

C^ ■ and cosmology. Among the various approaches to this problem (see for instance [l|), the use of the Dirac equation in 
an external gravitational field ([2|i [3|) has proven to be useful to obtain the semi classical corrections to the electron 
vj ' dynamics. However the results obtained in the litterature through a genuine Foldy-Wouthuysen transformation were 
^vj ■ restricted to a certain class of weak static gravitational field [3] and did not give any information about the dynamics 
,-_l ' in a general gravitational field. 

f^ , Moreover, these previous works did not take into account the role of Berry phases, which has recently proven to be 

\l ' very relevant in the context of the semi classical dynamics. Actually, since the seminal work of Berry [4], it is well 

^-:^ , known that a wave function acquires a geometric phase when a quantum mechanical system has an adiabatic evolution. 

(—K ' But it is only recently that a possible influence of the Berry phase on the transport properties (in particular on the 

— , ,. semiclassical dynamics) of several physical systems has been investigated. For instance, in [5| and [6[ the adiabatic 

'~i ■ evolution of the Dirac electron in an external potential was investigated. In [7|, the role of the Berry Phase on the 

' , Dynamics of an electron in a periodic potential has been explicitly derived in an Hamiltonian approach, in accord 

flN^" with the independent works of [8| and [9|. In our previous works, it was shown that position and momentum operators 

r^ , acquire spin-orbit contributions which turn out to be a Berry connections rendering the algebraic structure of the 

L| ' coordinates and momenta non-commutative (as found also in [9|). This drastically modifies the semiclassical equations 

of motion and implies a topological spin transport similar to the intrinsic spin Hall effect in semiconductor [10| . 

All these results presented similar patterns in the resulting coordinate algebra and equations of motion and it was 



X 



shown in llj that they all fall in a general formalism of semiclassical diagonalization for a large class of quantum 
system, including the Dirac Hamiltonian. This method has the advantage to derive directly the role of the Berry 
phase from the diagonalization procedure and introduces naturally the role of noncommutative dynamical variables. 

In this paper, we adapt the formalism developed in [111 ] to the case of a Dirac electron in an arbitrary static 
gravitational field. We find, at order ?i, the Foldy-Wouthuysen diagonalization for the Hamiltonian and assess the 
relevance of the Berry phases in the procedure. Our approach generalizes the results of [2| and [3|] (these last 
approaches, although different, being physically equivalent) since our method has also the advantage to be valid for 
all kinds of static gravitational fields and does not require a small field expansion. Moreover, we show that the 
Hamiltonian includes corrections similar to spin orbit and orbital momentum coupling found in [6] for Dirac electrons 
in electromagnetic fields. Here, the role of the magnetic field is played by the torsion of the gravitational field. 

We then derive the "covariant" coordinates algebra and their non commutative commutation relations. These 
variables differ from the ordinary position and momentum by the Berry terms and lead to different equation of 
motion than the usual ones [3|. We ultimately derive the semi classical equations of motion for the electron, showing 
the appearance at order % of new contributions to the force equation. 



ELECTRON IN A STATIC GRAVITATIONAL FIELD 



Consider an electron propagating in a static gravitation field where the corresponding metric satisfies goa — 0, 



which impHes ds^ = gooidx^)"^ — Oijdx^dx^. In the sequel, we will denote -^500 — l^(R-)j to be consistent with the 



notation of ('SJ], [2|)- Following \VM, the Dirac Hamiltonian for an electron in such a gravitational field has the form : 

H = y(R)a.P + ^e^^^rf (R)S^ + i^Tl^{IC)aP + 1/(R)/3m (1) 

and P given by Pc<=ft.^(R)(Pi + Tiegi3^^^{R)a'^) with /i^ the static orthonormal dreibein (a — 1, 2, 3), T^ ' the spin 
connection components and Ea/s-yCr'' = |(7"7^ — J^^")- It is known [12] that for a static gravitational field (which is 
the case considered here), the Hamiltonian H is herniitian. 

Because the components of P depend both on the operators P and R the diagonalization at order h is performed by 
adapting the method detailed in _lli to block diagonal Hamiltonians. To do so, we first write i/ in a symmetrical way 

in P and R. This is easily achieved using the Hermiticity of the Hamiltonian: H = ^ ( V^(R)a.P + P .aV{IV) 
jEp/^^Fn (R)S^ + F(R)/3to. The diagonalization is then performed in two steps, similarly to the method exposed in 

Diagonalization when P and R commute 

We consider a formal situation where R is first considered as a parameter r commuting with P. Some computations 
show that the Hamiltonian Hq (we add the index when R is a parameter) can then be diagonalized, at first order 
in h by the following unitary FW matrix 

C/o(P) = d(eo + V{r)m + ci/3 (v{r)a.P + P^.ay(r)) + Nj / ^J 2Eo {Eo + mV {r)) 



with E, = ^( v(^KP+f-^n.) ) + ^.^.(r), N = |H^I^, D = l+ |/3,|g^ and Fo, = e,,,Ff (r). 

The proof of this diagonalization relies on the fact that for each parameter r the matrices ft.^ and F" are independent 
of both the momentum and position operators, /3 and a.P anticommute and in the Taylor expansion of -Eo all 
terms commute with [3 and a.P + P+.a. In this context the diagonalized Hamiltonian is equal to UqHqU^ = 

^^( j(r)..P+P^.» VMy + ^2v2^r) +JtTo. (my(r)S + (sg;^^) with F,, = e,^,Ff (r), which reads 

U0H0U+ ^ (3^nVHr)g^Hr)P, + hs.p.rfh^^l^.P + m^VHr) + ^n. (mV{r)i: + ^J^^^l^^ (2) 

Corrections when P and R do not commute 

Now, we want to reintroduce the dependence in R, and apply the method given in [111 ] . Let us however remark that 
this method has to be adapted here since our Hamiltonian Eo(P) is not diagonal, but only block diagonal. Fortunately, 
given that the block diagonal part in Eo{P) is only of order Ti, our formalism still works at the semiclassical level. 

As a consequence, to perform our diagonalization procedure, it is sufficient to apply (at the first order in h) the 
following Foldy-Wouthuysen transformation 
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U(P,R)^ DlE + V{Ti)m+-f3[V{Il)a.P + F .aV(R)j + N j /.y2E {E + F(R)m) + X (3) 

and then project the transformed Hamiltonian on the positive energy states. The energy E in Eq|3] is given by 
expression Eq[2]where the parameter r is replaced by the operator R. All expressions in UiP) are implicitly assumed 
to be symmetrized in P and R. and the corrective term X must be added to restore the unitarity of U{P) which is 
destroyed by the symmetrization. |11| shows that this term is expressed as 

X = ^[Ap„Ar,\U{V,K) (4) 



where we have defined the position and momentum (non projected) Berry phases Ar — i?iUVpU~^ and Ap — 
—ihUy iiU~^ . The projection on the positive energy subspace that is needed to reahze the diagonaUzation leads us 
to define A^^ = P{Aii) and Ap = P{Ap). The resulting position and momentum operators can thus be written 
r — ihdp + A/j and p = P + Ap, where the explicit computation for the components of the Berry connections Ap 
and Ap gives 

^ 2E{E + my{Ii)) ^ ^ ' ^ ' 

^P- ^'' 2E{E + mV{IL)) ^''^^^ ^^> 

where E is the same as Eq above, but now R is an operator. Performing our diagonalization process ultimately leads 
us (after some computations similar to those presented in [l3|) to the following expression for the energy operator 
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e ~ £ + hB.'E/2e - (ARxp).B/e - ^^mhVV{r). (p x S) (7) 



where 




e = cJ{p, + h^^. { mV{r)i: + ,'"'''''', ,, ] ] VHr)g'^i (r) ( p, + h'^^. fmF(r)S + . ^^'^JP ,, ) ) + mW{r) 



{E + mV{r)) 



with Ti^ (r) — eafj-yT"^ (r) and the " magnetotorsion field" B is defined through the torsioni?-^ — —\PsT°'^^eai3^ 
where T^'P^ = /if (h^'^dih'''^ - h^f^dih'"^) + h^^vf^ - h'-'^Vf is the usual torsion for a static metric (where only space 
indices in the summations give non zero contributions) . Let us note that in Eqj5] we have neglected the curvature 
contributions that appear to be of order h . Interestingly, this semi-classical Hamiltonian presents formally the same 
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form as the one of a Dirac particle in a true external magnetic field [6|[ll|. The term B.u is the Stern-Gerlach effect, 
and the operator L = (^pxp) is the intrinsic angular momentum of semiclassical particles. The same contribution 
L appears also in the context of the semiclassical behavior of Bloch electrons (spinless) in an external magnetic field 
[7| y where it corresponds to a magnetization term. Because of this analogy and since T°'^^ is directly related to the 
torsion of space through T°'^^ = h^Ji^°'W^T^, we call B a magnetotorsion field. 

Let us also remark that the definition of the Berry curvatures (see below) , allows us to rewrite ultimately the energy 
e in a more compact form : 



^ = ^'\l[v^- ^.0-) V^v)g-^^ (r) [p, - ^^.@r^^ + mW{v)-^T,W^ (9) 

where 9'"'"'' = —^2ET~~r) ( 'TiT^(r)S-y + (E+mV(r)) 1 ^^ ^^^ "rescaled" Berry curvature. This formula clearly shows that 

is, the spin connection couples only to the Berry curvature. 

The semiclassical Hamiltonian Eql8] is the main result of this paper. It contains, in addition to the energy term 
e, new contributions due to the Berry connections. Indeed, the spin couples to the gravitational field through the 
magnetotorsion field B which is non-zero for a space with torsion. As a consequence, a hypothetical torsion of space 
may be revealed through the presence of this coupling. 

Let us note at this point that the coupling between torsion and spin has previously been studied by various authors 
([ij], [l5|, [16| for example and references therein), but in the context of general relativity and field equations, whereas 
the present paper consider gravitation as a fixed background. Our problematic is thus to describe the possible effect 
of torsion on the particle rather than deriving gravity field equations. From Eqs. [5]and[Bl we deduce the new (non 
canonical) commutations rules : 

[n,r,] = zhQlJ (10) 

[p,,p,] = iTief; (11) 

[Pi,rj] = -ihgij + ihO^^J (12) 
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where Qf = da^Ap^ — d^.Aa^ + [A^^.A/jJ. An explicit computation shows that 
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where the couple (P, R) has been replaced by (p, r) in P. We will need also 

-pjEi +p.S% 
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£(p,r)(£(p,r)+mV^(r))^'^'-'^~'' 
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Ultimately, using the commutation relations at hand plus the Hamiltonian, we can easily derive the equations of 
motion for the electron in a static gravitational field. As explained in ll[, the chosen dynamical variables are r and 
p. This choice comes quite naturally when considering the projection in the diagonalization process. It is also this 
choice that allows to take into account for the spin Hall effect for the photon, for example. We thus obtain 



r = (i-ep'')Vp£ + pxe'' 



-v^e.erf 



p = - (1 - e"') Vr£ + rxepp + Jvseef/ 



(20) 



However, these equations are incomplete per se since they involve the initial spin matrix ?iS that is not conserved 
through the dynamical evolution. To compute fully the dynamics for the electron at the first order in h, these 
equations have consequently to be completed with the dynamics of the spin matrix ?iS at the lowest order 



nt = l[s,.q=,^(2m(rar)xE) + 5(l4^)-PHpxi:; 
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The velocity equation contains in particular an anomalous velocity term pxG'"'', which has been described in several 
other circonstances. It is for instance responsible for the intrinsic spin Hall effect of Dirac electrons in electromagnetic 
fields |5|, 161]. In semiconductor, SO coupling being greatly enhanced with respect to the vacuum case, this anomalous 
velocity drastically modifies the transport properties of the charges [10| . A similar anomalous velocity appears also 
in the context of spinless electrons in magnetic Bloch bands [7] [8] . Here for the first time we predict an anomalous 
velocity contribution for an electron propagating in a static arbitrary graviational background. Note that this effect 
is present independently of the existence of a torsion of space. 

The force equation presents the dual expression f x 0pp of the anomalous velocity which is a kind of Lorentz force 
which being of order h does not influence the velocity equation at order h. It is interesting to remark that similar 
equations of motion with dual contributions pxO^^ and f x Q^^ were predicted for the wave-packets dynamics of 
spinless electrons in crystals subject to small perturbations [8]. 

Within our approach we can easily treat the ultrarelativistic limit iru? -^ which by simplifying all expressions 
gives a better feeling of the physics. 



Ultrarelatvistic limit 

In the ultrarelativistic limit rn — > 0, one recovers the same kind of Hamiltonian derived in |13| for the photon case. 
Indeed one readily obtain 

^^^^XpT^^X^B^ (21) 

where e — Ci (pi + A^IIlJ^ j g^Jg^^Q (pj + j jirj.p j ^-^^Yi X = p.Ti/p. Interestingly this energy can be expressed in 

terms of the helicity and not in term of E. In fact the terms j-V-Ee-Q^f and -^VseB^ in Eg EDI recombines with the 
gradient energy contributions which allows us to rewrite the equation of motion under the following form: 

f = (1 - 9j„.) Vpe + px9rr 

p = - (1 - Gp^) Vr£ + f X Opp (22) 

where in Vp£ we must consider p and S as independent variables. The term pxQrr causes an additional displacement 
of electrons of distinct helicity in opposite directions orthogonally to the ray. In comparison to the usual velocity 
f = Vp£ ^ c, the anomalous velocity term v^ is obviously small, its order w^ ~ cAV^jg*-' being proportional to the 
wave length A. 

To complete the dynamical description notice that at the leading order the helicity A is not changed by the unitary 
transformation which diagonalizes the Hamiltonian so that it can be written A = hp.Tl/p. After a short computation 
one can check that the helicity is always conserved 

for an arbitrary static gravitational field independently of the existence of a torsion of space. 

THE SYMMETRIC GRAVITATIONAL FIELD 

A typical example of such a metric is the Schwarzschild space-time in isotropic coordinates. This case, studied 
in a different manner in [2] and [3|. The Hamiltonian is given by, received a full independent treatment within our 
formalism in [11[. We thus present directly the needed results completed with the spin matrix dynamics. For a 
symmetric metric one has B.p == Fq = and the semiclassical Hamiltonian can take the following form [3|] 

i?o = ^ (a.PF(R) + F(R)a.P) + /3mV(R) (24) 

corresponding to the metric gij — 6ij ( -pTpr ) , <7io = and goo = ^^ (R-) • Similar computations to the ones performed 
in the previous section lead to the following expressions for the dynamical variables and the diagonalized Hamiltonian 

F2(R)S X P 
'2E{E + mV{R)) 
p = P+A+^P (26) 



r = R+yl+ = R-;i-_^A_A____ (25) 



e (p, r) ^ JF^(r)P^ + P'^F^M + mV^(r) - ^IJ mhW<j>(r). (P x S) (27) 

with (f> — j;. The commutators of the dynamical variables define the Berry curvatures [?'i,rj] = iTiQ^J, [Pi,rj] = 
—ihgij + ihQ^J and [Pi,Pj] = 0, with 
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2e3(P,r) V £(P,r)+my(r) 
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2,3(p,r)"^^'''^W(SxP)^. (29) 

(30) 



and 8[^ being unchanged and 8^ = 0. One can check, after developing r as a function of R and the Berry phase, 
that our Hamiltonian coincides with the one given in [3| when considering the semiclassical hmit (order h). This also 
confirms the validity of the Foldy Wouthuysen approach asserted in Q. The absence of dipole spin gravity interaction, 
in opposition with the transformation proposed in [2] , results only from a different choice of diagonalization matrix. 

Let us ultimately remark that, despite some similarities, our approach is different from the one developed in [3[. 
Actually, this last paper considers an approximated Foldy Wouthuysen transformation at the first order in a (weak) 
gravitational field. On the contrary, our approach, even if semiclassical, allows to consider a Foldy Wouthuysen 
transformation for an arbitrary gravitational field. 

To conclude this paragraph, we can easily derive the equations of motion from the commutation relations 

f = Vpe-Pxe'^'^ + JVse.Olf (31) 

P = -Vr£ + VrE.eP'' (32) 

= _Vre + Vr,£^^^^mV.0(r) (P x S)^. (33) 

nt = "^^W^ sx (VV^(r)xP) - -Ex (Vi^(r)xP) 
e{e + mV{r)) ' ' ' ' e 

Although the first equation differs from the one ontained in [3j due to a different choice of dynamical position variable, 
the last two equations reduce to [3], in the case of weak fields, : 

p = _!!^vnr)-^Vi^(r) + -/^V(S.(Vnr)xP))-Av(S.(Vi^(r)xP)) 

' TTi n f) 

n-E = ^Sx(VF(r)xP)--Ex(Vi^(r)xP) 



where here e = \JjP^-Vrr? . 

Let us stress again that our dynamical position variables are not the canonical variables R chosen in [3| , but the non 
commutative ones r. As a consequence, our equations obviously differ from theirs by the role of the Berry curvature, 
even if the Hamiltonians are formally identical. 

CONCLUSION 

The semiclassical limit for Dirac particles interacting with a static gravitational field was investigated. We found 
a Foldy- Wouthuysen transformation which diagonalizes at the semiclassical order the Dirac equation for an arbitrary 
static space-time metric (not necessarily symmetrical contrary to previous works). The main results are the following. 

First, the energy includes, through the Berry phases, a coupling between the torsion of the space and the spin of 
the particle. Second, we have been able to derive the semi classical equations of motion for the electron. Our set 
of dynamical variables is non commutative and this choice induces new corrections to the dynamical equations that 
could be relevant in a strong gravitational field. Ultimately, we retrieve the particular case of a symmetric static 
spacetime metric [3| , and confirms the generality of our approach 
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